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The ground-state properties of neutron-deficient Hg isotopes have been investigated 
by the constrained self-consistent Hartree-Fock-Bogoliubov approach with the Skyrme 
and Gogny effective forces. In the case of the Skyrme interaction we have also applied 
the Hartree-Fock+BCS model with the state-dependent 5-pairing interaction. Potential 
energy surfaces and pairing properties have been compared for the both types of forces. 



1. Introduction 

The self-consistent mean-field Hartree-Fock-Bogoliubov (HFB) method is one of 
the standard approaches in nuclear structure and low-energy dynamics research.^ 
In this paper, we compare two self-consistent models: the HFB framework with 
the Skyrme interaction^ and the HFB theory with the Gogny forced We have 
tested these two models with various parameter sets. For calculations we have 
chosen the commonly used SkM* and SLy4 parameterizations of the Skyrme force. 
We have also performed calculations with the SkM* parameterization the Hartree- 
Fock+BCS (HF+BCS) model. This has been done in order to check the impact of 
a non-self-consistent BCS scheme on the results in the case of nuclei far from the 
stability line. The acknowledged D1S parameter set of the Gogny force has been 
set together with recently developed DIN and DIM sets. 

The neutron-deficient nuclei, in the vicinity of the neutron mid-shell at iV=104 
and near the Z=82 shell closure, exhibit a variety of coexisting shapes (see, e.g., 
Refs.EEEl). 

The potential energy surfaces (PESs) of nuclei from this region show 
complex structures around their ground states. This brings about a stringent re- 
striction on the validity of various nuclear models and forces parameters. We have 
chosen light 178 . 18 °. 182 . 184 Hg isotopes for tests of the mean-field models with the 
Skyrme and Gogny forces. 
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2. Calculation details 

Calculations with the Skyrme interaction have been carried out using a symmetry- 
unrestricted code HFODD (v.2.38j jZISIHI, that solves the self-consistent HF+BCS 
and HFB equations in the Cartesian deformed harmonic-oscillator basis. In the 
calculations, for the basis we have taken the lowest 1140 stretched single-particle 
states originated from the lowest 31 and 26 major oscillator shells for the HF+BCS 
and HFB approach, respectively. 

In the p/i- chan nel we have applied the Skyrme energy density functional with 
the SkM* RefP^ and SLy4 RefP^ parameterizations. In the pp-channel we have 
used the MIX variant of state-dependent (5-interaction (see, e.g., Reff™ 
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where po = 0.16 fm and a = 1. 

The pairing strengths for neutrons/protons Vq" — —372.0 MeV, 
-438.0 MeV (SkM*-BCS), 7 " = -268.9 MeV, V P = -332.5 MeV (SkM*-HFB), 
and Vo" = -308.0 MeV, V P = -343.4 MeV (SLy4-HFB) have been adjusted to 
reproduce the experimental pairing gaps in 252 Fm. As we deal with contact in- 
teractions, for the HF+BCS approach we have used a finite pairing-active space 
defined by including N/Z lowest single-particle states for neutrons/protons, while 
for summing up contributions of the HFB quasi-particle states to density matrices 
the cutoff energy of 60 MeV has been used. 



Table 1. Gogny parameters in D1S, DIN and DIM parame- 
ters sets 



parameter 


D1S 


DIN 


DIM 


Wx [MeV] 


-1 720.30 


-2 047.61 


-12 797.57 


Si [MeV] 


1 300.00 


1 700.00 


14 048.85 


#1 [MeV] 


-1 813.53 


-2 414.93 


-15 144.43 


A/i [MeV] 


1 397.60 


1 519.35 


11 963.89 


Ml [fm] 


0.7 


0.8 


0.5 


W 2 [MeV] 


103.639 


293.02 


490.95 


B 2 [MeV] 


-163.483 


-300.78 


-752.27 


H 2 [MeV] 


162.812 


414.59 


675.12 


M 2 [MeV] 


-223.934 


-316.84 


-693.57 


M2 [fm] 


1.2 


1.2 


1.0 


i [MeV fm 3 ( 1+ T)] 


1 390.6 


1 609.46 


1 562.22 


x 


1 


1 


1 


7 


1/3 


1/3 


1/3 


W LS [MeY fm 5 ] 


130 


115.0 


115.36 



The most popular parameter set of the Gogny force D1S was fitted almost 
twenty years ago ff^l Recently, two new parameter sets have been proposed: DHsfl^ 
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and D1M.Q23 These parameterizations have been developed with special care about 

accurate reproduction of nuclear matter properties and dependence of symmetry 

energy on nuclear density. The DIM set has been fitted to reproduce over 2000 

experimental masses after including beyond-mean-field quadrupole corrections. All 

parameters of the conside red Gogny forces are given in Table [5] 

FTTH 

The Gogny calculations u have been performed in the deformed axial basis with 
N z — 22 oscillator shells in z direction and N± — 15 in perpendicular directions. 
Oscillator lengths have been optimized in every point of calculations. 



3. Results 
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Fig. 1. (Color online) Potential energy curves from the Skyrme and Gogny self-consistent mod- 
els restricted to axially symmetric shapes as a function of Q2 mass moment calculated for 
178, 180, 182, 184 Hg i aotopes (f or details, see text). 
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In Fig. [T] we present potential energy curves of four neutron-deficient isotopes 
I78,i80,i82,i84jjg ca i cma t ec i f or axially symmetric shapes for moderate deformations 
(IQ2I < 50 b). The energies have been calculated with the use of two different ef- 
fective forces: Skyrme and Gogny. In the case of the HFB-Gogny approach we have 
used the DS1, DIN and DIM parameter sets. The SkM* and SLy4 Skyrme param- 
eterizations have been applied to the HFB-Skyrme framework (curves marked as 
SkM*-HFB and SLy4-HFB in Fig.HJ). In addition, the SkM* set has been applied to 
the HF+BCS model with the MIX paring S- interaction (SkM*-MIX). The potential 
energy profiles shown in Fig. [T] include only the mean- field energies with projection 
onto good center-of-mass momentum. Other beyond mean- field corrections, e.g., a 
rotational energy correction, have been neglected. 



(SkM* - MIX) (SkM* -MIX) 




OOOOOO OOOOOO 

^ ki ui !fc ui j k) w i in bi 



P P 

Fig. 2. (Color online) Relative potential energy surfaces calculated in the SkM*-MIX model for 
178,180,182,184 Hg in £ _ 7 p l ane . 



In spite of the fact that we have used various mean-field approaches, results 
obtained with all parameter sets are quite similar. First of all, one can easily find 
five minima in all presented isotopes and models: two at the oblate side (Q% « 
—20 and —10 b) and three with prolate deformation (Q2 ~ 5, 20 and 40 b). The 
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most pronounced prolate minimum is localized at Qi ~ 20 b and the HFB models 
predict ground state there. Nevertheless, in all considered nuclei, the deepest oblate 
minimum at Q2 ~ — 10 b is not more than 1.2 MeV higher for the Gogny forces, 
and even less for the Skyrme forces. It is worth to note that the deformed prolate 
character of the ground state is best seen in 180 Hg and 18 2 Hg for all models. This is in 
accordance with experimental results published recently--. In 178 Hg and 184 Hg the 
HFB models predict oblate minimum on 0.5 MeV above the prolate ground state, 
whereas the SkM*-MIX model indicate oblate ground state. The other minima 
(Q2 w ^20, 5 and 40 b) are not distinct, especially in isotopes 182 Hg and 184 Hg 
where they convert into shoulders only. 
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Fig. 3. (Color online) Neutron and proton pairing energy for isotopes and models presented in 
Fig. 1. 
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The PESs calculated with the DIN and DIM Gogny forces lie higher by about 
6 MeV than the values obtained with the D1S force. The binding energies in 
these two parameterizations coincide for spherical shapes and do not differ more 
than 1 MeV for small quadrupole deformations. Inclusion of the beyond-mean- 
field quadrupole corrections in the DIM set improve substantially agreement be- 
tween theoretical and experimental masses ff^ However, in the case of four neutron- 
deficient Hg isotopes investigated here, the D1S parameter set gives the binding 
energies closest to the experimental values, marked in each panel by horizontal 
lines. 

The PESs calculated with the Skyrme SLy4-HFB model are shifted up by around 
3 MeV in relation to the values obtained with the D1S Gogny force. In both the 
SkM*-HFB and SkM*-MIX models the calculated PESs are almost identical, how- 
ever, their binding energies overestimate the experimental values by around 10 MeV. 
This is the worst agreement for all models investigated here, which indicates that 
the SkM* set fails in prediction of absolute experimental data for nuclei in this 
region. 

For better understanding the relations between minima presented in Fig. 1, we 
have performed calculations with broken axial symmetry. Fig. 2 shows the relative 
potential energy surfaces calculated in the SkM*-MIX model for 178 . 180 > 182 > 184 Hg in 
the (3 — 7 plane. For 180 > 182 Hg isotopes, one can find shallow prolate global minima 
at f3 « 0.3 that correspond to the ground states at Q2 — 20 b in Fig. 1. However, 
in the case of 178 < 184 Hg isotopes, the SkM*-MIX model suggests even more shallow 
oblate ground state minima at (3 « 0.15 (Q2 = —10 b). In the non-axial calculations 
these oblate minima have appeared to be extremely soft toward the 7— deformation 
and one can find no barrier that separates them from the prolate minima at (3 « 0.3. 

In order to give further insight into the relation between the mean-field models 
we have provided neutron and proton pairing energies of Hg isotopes in Fig. 3. The 
results presented here have been calculated in the same frameworks as in Fig. 1. 
One can see that almost all results coincide, only absolute value of proton pairing 
energy in the SkM*-HFB model is larger. 

4. Conclusions 

We found, as expected, that various mean-field parameterizations give slightly dif- 
ferent absolute values of the total energy. However, it is encouraging that they 
give quite similar dependence of the energy on the deformation. The results of 
this research indicate also that despite of differences between considered mean-field 
models the predicted pairing properties of investigated light Hg isotopes are very 
similar. 

An important issue, not addressed in this study, is a comparison with experiment 
which needs a thorough consideration of various correction terms such as those 
resulting from the particle number or angular momentum projection. The results 
reported in Section 3 imply that it is necessary to include non-axial deformations 
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in the case of considered Hg nuclei. Because of large 7— softness of the potential 
energy it does not seem sufficient to consider only pure prolate-oblate shapes. 

Finally, recent achievements in the self-consistent treatment of quadrupole col- 
lective excitations (in the spirit of the Generator Coordinate Method or the Adia- 
batic Time Dependent HFB method) open a new possibility to test the mean-field 
models using rich experimental data on low-lying collective states. 
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